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Abstract 



In this paper we investigate the smallest eigenvalue, denoted as Atv, of a (iV+l) x (JV + 1) 
Hankel or moments matrix, associated with the weight, w(x) = exp(— x@), x > 0, (3 > 0, in 



the large N limit. Using a previous result, the asymptotics for the polynomials, P n (z), z ^ 



[0,oo), orthonormal with respect to w, which are required in the determination of \n are 
found. Adopting an argument of Szego the asymptotic behaviour of Atv, for /3 > 1/2 where 



the related moment problem is determinate, is derived. This generalises the result given by 
Szego for = 1. It is shown that for (3 > 1/2 the smallest eigenvalue of the infinite Hankel 
matrix is zero, while for < [3 < 1/2 it is greater then a positive constant. This shows a 
phase transition in the corresponding Hermitian random matrix model as the parameter f3 
varies with (3 = 1/2 identified as the critical point. The smallest eigenvalue at this point is 
conjectured. 



1. Introduction. 

In the theory of Hermitian random matrices, the Hankel determinant plays an important role, 



Dn = oMn^ ■ 



For a given weight function w(t) on J (C _R,) the moments fi^ are 
1.1) /i fc := / w(t)t k dt ; A; = 0,1,2,- •• 



Associated with w(t) is a Hankel matrix or moment matrix of order N + 1, {Hjk}, whose entries 
are given by 

(1.2) H jk := N+k ■ 0<j,k<N. 



It is believed that correlations between eigenvalues of random matrices are universal after a 
suitable rescaling. In the following treatment we will show that a fundamental quantity, namely 
the least eigenvalues of these Hankel matrices exhibit a critical dependence on the weight function. 
It is this non-universal property that motivates our investigation of this problem. 

If J is a single interval say [a, b], where a and b are fixed and the Szego condition, 



v (x)dx 



< oo, v :— — \nw, 



(b — x)(x — a) 



is satisfied then the asymptotic behaviour of the Hankel determinants for large N was established 



by Szego, fllP] . Let Xn denote the smallest eigenvalue. Szego also investigated the behaviour of 
Xn for large iV ||. He studied the cases for which J can either be a finite or infinite interval 
with special choices for w. If w(x) — 1, x G (—1, 1) and w(x) = 1, x e (0, 1), then the respective 
smallest eigenvalues are for large N [] 



\2N+3 



X N ~ 2^N^(V2-1) 

X N ~ 2 i MiV^(v / 2-l) 4JV+4 . 



Widom and Wilf [TT| generalised Szego's results to a kind of "universal" law. Thus if w(x) > 



0, x G [a, b] and the Szego condition is satisfied then it was found in |]TT[ that 



A 



N 



AN^ B 



-N 



where A and B are computable constants depending on w, a, b and are independent of N. 

In 0, Szego also considered the cases of infinite intervals where w(x) = exp [— x 2 ] , x e 
(— oo, +oo) and w(x) = exp [—x] , x G [0, +<x>), are the weights of the Hermite and Laguerre 
polynomials^. The respective smallest eigenvalues are 



Xn 

Xn 



2-71-3 eA^exp -2(2N)^ 



22 7T 2 e iV^exp 



-47V2 



1 Throughout this paper, the relation, ajy — bjy means limjv->oo «» /6jv = 1- 

2 There is a factor of 4 missing from the original formula for Ajy; the last equation in page 677 of 



Observe that in the examples given above the smallest eigenvalues are exponentially small. There- 
fore it is very hard to numerically invert the Hankel matrices associated with these weights. 

It is well known that A at is given by the Rayleigh quotient: 

(1.3) Ajv = muW ■>' \ . 

{ l^j=0\ x j\ ) 

If ttn(z) is a polynomial of degree N, with coefficients Xj, j = 0, ..., N 

N 

(1.4) *•*(*):= 5>,V , 

3=0 



N - f 

1.5) Yl H jk XjX k = / \-K N (t)\ 2 w(t)dt 

■ 1. n J J 



then 

N 

E 

j,k=0 

and 

(1-6) El^| 2 =rMe^)|f 

Consequently we can rephrase the extremal expression for \ N , (1.3), as, 

2tt r r 2 ^ 



1.7) ^l = max <{/ |7TAr(e^)| 2 # : / |7TAr(t)| 2 «;(t)dt = 1 
Aat Uo Jj 

Let {P ra (t)} be the polynomials, orthonormal with respect to w (£), then -jin has the expansion, 

N 

(1.8) 7T^)=E C i P i(«)- 

j=0 



Thus 

E 

i,fc=0 

where 

r27r 



/■27T W 

;i.9) / |7r w (e^)| 2 rf0= X! K ]k cfc k 



1.10) K jfe := / P^P^z)^ ; z = e l 

J U 



Therefore (1.7) is equivalent to 



N N 



;i.ll) — = max<{ Y, K JkCjC k : E l c i 

J,k=0 j=0 



Xn 



With the Schwarz inequality, which states that for all values of j and k 

\K jk \ < K%K\ k , 
and Cauchy's inequality we obtain an upper bound of (1.11): 

TV N 

(1.12) Yl K jk CjC k < Y \ K jk\\cj\\ck\ 

j,k=0 j,k=0 

N I I 

< Y KlKl k \c,\\c k \ 

j,k=0 

< (e^) (l> 

\j=o ) \j=o 

N 

= E^i- 

Therefore a lower bound for the smallest eigenvalue Xn is given by 

(1-13) ^j-<X N . 

2^j=0 ^33 

This paper is organised as follows: In section 2, by adopting a previous result 0, we obtain the 
asymptotic formula for the polynomials orthonormal with respect to w(t) := exp[— t^], f3 > 1/2, 
which is then employed in sections 3 and 4 for the determination of the large N behaviour of 
Xn- In these sections we show, following 0, by an appropriate choice of the vector {cj}, that the 
lower bound given by (1.13) is in fact an asymptotic estimate for large N. By a simple application 
of Laplace's method, J2f=o Kjj is estimated. Thus the asymptotic form of Xn follows. In order 
to test the accuracy of the theory, these results are checked against numerical calculations for 
various /3 and N, which were obtained using the Jacobi rotation algorithm [Q to reduce the 
Hankel matrix to diagonal form. This is found in section 5. 

2. The weight w(t) = exp[-tP], t e [0, oo). 

In this case, the moments are 

(2.1) *. = f[—)- 



In order to find a lower bound for the smallest eigenvalue good knowledge is required of the 
associated orthonormal polynomials {P N (z)}, for N large and z £ (0,oo). In ||, by applying 
the linear statistics formula for matrix ensembles together with the Heine's determinant repre- 
sentation, asymptotic forms for the polynomials with weight w(t) = exp[—v(t)], where v(t) is an 
arbitrary convex function supported on [0, oo), are derived. The zeros of these polynomials are 
supported on (a, b) C R. Here a = 0, whilst b(N) follows from the condition that ensures that 
P/v(£) has N roots on (a, b), one finds that ||, 



(2.2) b(N; p) = CN^, where C = C{(3) : = 4 



r 2 (/3) 



i 



r(2/?)_ 

The normalised polynomials as N — > oo are found, using ||, to be 

K J w V2ri [C(i + C)]* b 

where / is given by 

(2.4) flO.&W , ( ^,. 

2tt Jo y-tjyty-y) 

From the definition and basic properties of the hypergeometric functions [|7| , 

(2.5) f{t) = -^TV / C( 1 + C) 2 Fi(l,l-/3^-/5;-c)-^sec7r/5 



N C f 1 n 1 

2 F 1 [ !,-;/?+ 1; 



/3 V l + C V 2 1 + C/ 

At this point note the dichotomy of the problem, the nature of the Hypergeometric function 
dictates that whilst the first representation is more convenient in the large b limit, where |£| << 1, 
it cannot be used when (3 — n + ^, n = 1, 2, . . ., necessitating the use of the second result of 
(2.5) in such instances. 

Using the fact that 

1 1 3 



2-6) HfC + Jl + C) = fC2F 1 (-,-;-;-C) 



we find, 



t) 



p E \P~h\ 



2.7) (2iV + l)ln(^CH-V 1 + C)^V?^ £ (-^rC 

V 7 ™-^ r =Q 

5 



\ r n n'+^-i 3 



2 



where E[n] denotes the integer part of n and 



(2.8) a r :-- 



T(r 



So the asymptotic expression of the polynomials for t <£ (0, oo), is, 



''-1V V < T / (-f)l 3 E[P ~^ ] 

(2.9) P N (t)^^d=k r exp\-f(t) + ^ 7 = £ (-l) r a r C+i-'- ! 



V=27rt 



C?V* £> 



To make further progress we now consider separately the two possible cases, as identified above, 
for (3 > 1/2. 

3. (3^n + \, n = l,2,3--- 

When /3 ^ n + |, we use the first form for f(t) in equation (2.5). The series expansion for 
the function 2 Fi(i,i-/3;|-/3 ; -c), valid for \(\ < 1, is 

r(|-/g) f r 1v . r(i-/? + r) 

T(l-/3)^ j r(|-/3 + r)^ ' 



(3.1) 2 F 1 (l,l-/3;--/3;-C 



whilst for |£| < 1, -y/1 + C may be written as 
i _i °° Y(r - -) 

With this noted, the expansion for f(t) as £ — > is 



(3.3) /(;) 






2^r(i- i) VcJ S ( iy7vC 



:-*y 



sec7r/5 



where 

(34) t . = ^ r(.-i)r(i-j? + r-.) 
lJ ' 4 ' 6 " ir„ r(» + i)r(| - ,3 + r - S ) 



Recall that £ = — tC X N p , and by the use of equation (2.9) we have, 



(3.5) P N (t) 



\N 



^7T 



"*)' 



sec7T/9 



x exp 



(-tCN^)~exp 



with 
(3.6) 



l\ T '. Chrp ~~T 



2T(l-f3) 



Note with (3 — 1, we find C = 4 and Aq = 4y/n and consequently recover the classical result 
for the Laguerre polynomials due to Perron f9j, 



(3.7) P N (t) 



-1 



>iv 



2A 



{—tN) * exp 



2V-tN + 



t<£[0,oo), N ^ oo . 



With P/v(£) having the form (3.5), where A = -£0§ is positive for (3 > 1/2, we observe that 
for sufficiently large j and k the dominant contributions to Kj k are from the arc of the unit circle 
around t — — 1. Thus by fixing an arbitrary positive number u and confining ourselves to values 
of j and k satisfying 

(3.8) N - uNWi <j,k<N, 
we have 

(3.9) K jk ~ jT^ P, (e^) P fc (er*) d</> . 

Using the substitution 9 = <ft — tc and expanding the integrand for \9\ << 1 gives the following, 

(3.10) K jk : 



(-l) j+k i 

2ttVC ' 



X 



cxp 



1 -if 2 j 4 

_^ y (-i) r — 

(2r + l)i9 , . x _± 



2n2" 



(2r + 1) 2 9 



(I 1 _r -j 1 r_\ 

j 2 ' 3 9 + k ^ J 



(i 1 r_ i 1 _r\ 

j 2/3 P - k 2,3 ,dj 



d9 . 



1 1 i_r -I 1 r_ 

Because j 2 ' 3 p — k 2 f> p remains bounded in the range specified by (3.8) we can disregard 
the linear term in 9 in the integrand. This integral can then be approximated by extending 
the range of integration to the real axis, which does not affect the asymptotic behaviour, as 
contributions from (— oo, — e) and (e, oo) are sub-dominant compared to those from [— e, e] as 
j, k — *■ oo. Therefore, 

(3.11) K jk ~ \ ^j AoV^iV^-ig 



(ttC)4 
x exp 



1 B P-3] A 

= E (-i) r £ (/"*-* + * 1- ' 



c 



r=0 



J r_ 

2/3 ^ 



From (3.11), we see that when j and k are sufficiently large and satisfy (3.8), 

(3.12) K jk ~(-lY+ k K} J KJ k . 

This is especially useful as it enables the determination of the large N behaviour of Ajv- By 
choosing the vector {cj}, as in ||, such that 



(3.13) Cj 



ae inj K^ if E[N - uN^} <j<N 



if j < E[N-uNw] , 

where a is a positive number determined by the condition 



N 



(3-14) £ 



Of 



3=0 



N 

° 2 E 

i 

j=B[iV-a;iV2?] 



*ii = i 



we find, using (3.12) and (3.14), that 



N 



N 



(3.15) Yl K jkCjC k 

j,k=0 



1 1 

2 

A7v 



i 

j,fc=£[7V-w7V^] 



<7 



# 



TV 

E 

i 

j=E[JV-wJV2?] 

V 



.7.7 



E 

i 



Recalling equation (1.11), we see that since u is arbitrarily large the asymptotic behaviour of 
the maximum, by virtue of the inequality (1.13), is well approximated by ^2j =0 Kjj. Therefore 
we have shown that 



2ty 



N 



(3.16) —^J2 K n- 

The leading behaviour of this sum for large N is in turn found by replacing the sum by an 
integral and by applying Laplace's method, which in this context may be stated as : 

If for x G [a, b], the real continuous function <f>(x) has as its maximum the value <f)(b), then as 

N ^ oo 



(3.17) / f{x)exp[N<f){x)}dx 

J a 



/(6)exp[iV0(6)] 
N<t>'(b) 



A simple calculation gives the expression for A 



N, 



(3.18) -^ ~ -n~^C^Aje secwl3 N-^ + ^ 
An 4 



H) 



x exp 



2N l ~^ E[ ^ ] A 

Z JL = . y (_i)^iv-f 



r=0 



Putting [3=1, Szego's classical result for the Laguerre weight is recovered: 

(3.19) -^ ~ 2 _ t 7r -| e - 1 Ar-3 exp^V^V] . 
an 

From (3.18) we see that the smallest eigenvalue is exponentially small for large N and is zero for 
the corresponding infinite Hankel matrix. 

4. p = n+l, n = 1,2, •■■ 

In this section we investigate the case where j3 = n+ \ , n > 1. Such cases, as was explained 
previously, require the second form of fit) in (2.5). To obtain the asymptotic expansion for fit), 
we first note the following result for the hypergeometric function : 

If /3 = n + \ with n = 1, 2, . . . then 



(4.1) 2 F 1 [l,-;(3 + l;x 



l)f>-$ 



13+1 



v^rln 



X' ' 2 



l±Vx 
1-v/i 



+ 



M 



r =i L r _i \x - I 



where L r is given by 

(4.2) L r := ^-<T(r) . 

This is easily be proved by using an inductive argument, noting the following version of Gauss' 
recursion relations Rl 



!!■•>) ,/•', (l,I;n+|z 



+ 



(n + §)(*-!) 
(n + l)z 

n(n + |) 



1 3 

2 Fi [l,-;n + -;z 



(n+l)(n+i) 



2' 2 

1 3 

2' 



2 Fi (1, -;n + -;^ 



Y- 2-Fi ( 1, -; ri + -;z 



"'2' 



together with the fact that 



(i- 1) ^W 1 ^!* 



3(s-l) 



In 



z-i 



1 + A/i 



1-x/i 



3 



Therefore, 



(4-5) f(t) 



-1) 



p+i 



2n 



-tf In 



v / TTC + i 
v / T+C-l 



l + C 



E( 

r=l 



rC" 



L„. 



Using (3.2), we find 



(—■\\P+h(—t\P 
(4.6) f(t) ~ l ; _ l ; hi 



-t) 



2ti 



4tt2 r=0 



r+i-/3 



ICI«i 



where 



-M 



7« 



(4-7) 5 r := £ 

8 = 1 L S-I 



and 



f r(r-i/2) if r > n 

(4.8) 7r := ^ r(r+i) "^ 

I if r < 0. 

Recalling £ = —tC~ l N~~P, the strong asymptotics of the polynomials for £ ^ [0, oo) reads, 

(4.9) P N (t) ~ 



^A ? 



J) -tCJV3)*exp 



'27T 



(_l)/3-i(_ t )/3 / 4C <]v? 
— In 



2^ 



-t 



x exp 



N 1 '^ Ij 



vc 



E(-i) r ^ 



:-*) 



r+J 



(CW 



o V 



where 
(4.10) B r :-. 



2P 6 ^' 



Note the appearance of the logarithm in exponential. Since -Bo = 2 g~i > ® anc ^ usm g an 
argument similar to that in the previous section, we see that in determining Kjk the essential 
contribution comes from the arc in the vicinity of t — — 1. As before restricting j, k to the range 
given in (3.8), we have, 



(4.H) K. 



jk 



Pj{-e^)P k {-e-^)dB . 



10 



We expand the exponential in the integrand for \9\ « 1, keeping terms up to second order and 
then extend the range of integration to the infinite interval 
ln(j/k) remain bounded in the range given by (3.8), we find 



-| _ 1 r_ -I 1 r_ 

then extend the range of integration to the infinite interval. Because j w p — k 2 (> p and 



(4.12) Kj 



;M : 



{-i) j+k -i ii, i, (-if-? 



(ttC) 
x exp 






Again note that for sufficiently large j and k, satisfying (3.8), 



(4.13) i^ (-1)^44 



l ij- tv fcfc 



Repeating the argument of the previous section, it follows that 



(4.14) — ~ / K 31 dj . 



The leading term in the asymptotic expansion of this integral as N 
application of Laplace's method and is given by 



oo follows from an 



(4.15) 



2vr 
Xn 



1 li-i 1,1, i (-i)' 3- ^ 



x exp 



9/V 1- ^ ^ 2 Br 



r=0 



(-1)0-1/2 

Effectively exp [sec it (3] in (3.18) is replaced by (ACN 1 ^) * . Note the alternating nature 
of this additional factor depending on whether (5 — 1/2 is odd or even. Again (4.15) shows 
that limAT^oo Aat = 0. According to standard theory [[j]], the moment problem associated with 
w(x), x > is indeterminate if 



v(x) 



\/x(l + x) 



dx < oo. 



Therefore j3 = 1/2 is special as it marks the transition point at which the moment problem 
becomes indeterminate. Assuming, the result given in (2.9) holds, we have 



(4.16) P N (t) 



\N 



2tt 



r —t) *N 2 exp 



In 



7T 



AirN 



+ 1 



t$ [0,oo). 



11 



Again if we confine ourselves to the range where j and k are sufficiently large to enable the use 
of the above asymptotic representation, we find that the major contributions to K 3 ^ are from 
the arc around t = — 1. But, due to the behaviour of P/v(t) with increasing N, it is quite clear 
that \Kjk\ decreases as j, k — > oo, making an analysis analogous to that of the previous sections 
impossible. 

It is however possible to obtain an approximate lower bound for the least eigenvalue, since 
(1.13) still holds. Applying the Christoffel-Darboux formula |J and the result given in Q for 
the large N off diagonal recurrence coeeficients, we find, 

N „„ N 



(4.17) Y. K H = I H p A-e i6 )PA-e- i6 )de 

„'_n J— IT _•_ n 



j=Q ->-*j=0 



r^j 



IX 



-2 at2 



- P N (-e w )P N+1 (-e - ie ) - P N { -e^ e )P N+1 {-e^ 

e i9 _ e -i& 



N 2 / JVV y JV ^ 1V y JV ^ 1V J d6 



Thus using Laplaces method, 



b / 27T 

dx f(x)exp[N(j)(x)] ~ f(c)exp[N(f)(c)] J _ N , as A^ 

where c G (a, b) is the maximum of <f>{x) for x G (a, 6), gives 



(4.18) E^. 



n 



(AixNef 



j=o A^J\n(ATcNe) 

So at the point j3 = 1/2 the smallest eigenvalue appears to decrease algebraically instead of 
exponentially. 

5. Numerical Results 

In this section we check the accuracy of our asymptotic expressions for the least eigenvalue 
of the the various Hankel matrices against numerical results. Due to the fact that the moment 
matrices in these cases are very ill conditioned becuase of the vast range in scale of the matrix 
elements, the Jacobi rotation algorithm [|i~2"] , proved far more stable than the more conventional 
techniques for numerically determining a small selection of the eigenvalues of large symmetric 
matrices such as the Lanczos procedure or Householder's method ||. This appears to be an 
unusual phenomenon. Because of the behaviour of the matrix elements in these problems it 
is necessary to implement a multiple-precision package that allows floating point arithmetic of 
arbitrary precision. The library of sub- routines created by Brent was employed to combat the 
effect of rounding errors in the numerical procedures. 



12 



For < (3 < 1/2, the corresponding moment problem becomes indeterminate [l], and as a 
consequence the sum 

OO 

£l^(*)f> 

converges for every z in every compact subset of the complex plane. Therefore 

oo 

and the smallest eigenvalue for the corresponding infinite Hankel is a positive constant bounded 
below by 27r/£. Proof of the extention of the above statement to all indeterminate moment 
problems and other related topics can be found in g. The situation for < (3 < 1/2 is in 
contrast to the results for (3 > 1/2 where (3.18) and (4.15), as confirmed by the numerics, show 
that the sum diverges - A fact that is also well-known from the standard theory when the moment 
problem is determinate M, This separation of behaviour in the two regions is the phenomenon 
of phase transition alluded to earlier. 

The comparison between the numerical values of A n and those obtained from the theoretical 
expressions (3.18) and (4.15) is shown in table 1. and figure 1. below. 



300 



Figure 1: The percentage error of the theoretical values of A at when compared to those obtained 
numerically, for various (3. 
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p 


iV 


Numerical Atv 


Theoretical A^r 


1 


50 


2.0948 x 10" 10 


2.3695 x 10" 10 




100 


2.1079 x 10" 15 


2.3006 x 10" 15 




150 


2.9551 x 10" 19 


3.1743 x 10" 19 




200 


1.6387 x 10~ 22 


1.7437 x 10~ 22 




300 


5.5215 x 10" 28 


5.8090 x 10" 28 


3 

2 


50 


6.4066 x 10~ 22 


6.8438 x 10~ 22 




100 


6.2353 x 10~ 36 


6.5384 x 10~ 36 




150 


9.9476 x 10~ 48 


1.0343 x 10" 47 




200 


2.8132 x 10~ 58 


2.9101 x 10" 58 




300 


4.6009 x 10" 77 


4.7300 x 10" 77 


7 
4 


50 


6.4483 x 10" 27 


6.6844 x 10" 27 




100 


1.6976 x 10" 45 


1.7424 x 10" 45 




150 


1.5193 x 10" 61 


1.5525 x 10" 61 




200 


3.9265 x 10" 76 


4.0009 x 10" 76 




300 


1.4844 x 10~ 102 


1.5074 x 10~ 102 


2 


50 


2.7356 x 10" 31 


2.5449 x 10" 31 




100 


3.8907 x 10" 54 


3.6415 x 10" 54 




150 


2.9557 x 10" 74 


2.7769 x 10" 74 




200 


8.9775 x 10~ 93 


8.4574 x 10~ 93 




300 


9.5593 x 10~ 127 


9.0396 x 10~ 127 


5 
2 


50 


2.2384 x 10~ 38 


2.4010 x 10~ 38 




100 


1.2580 x 10" 68 


1.3288 x 10~ 68 




150 


5.3195 x 10~ 96 


5.5789 x 10~ 96 




200 


1.2155 x 10~ 121 


1.2691 x 10~ 121 




300 


1.5236 x 10~ 169 


1.5819 x 10~ 169 



Table 1: Numerical and theoretical values of A^ for various /3 
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